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^ ■ Abstract 

We consider the downlink of cellular systems in which the base station employs multiple transmit 
antennas and each user has one receive antenna. We consider communication schemes in which the users 
have Quality of Service (QoS) requirements, and we study the design of robust broadcasting schemes that 
minimize the transmission power necessary to guarantee the QoS requirements for all channels within 
bounded uncertainty regions around the transmitter's estimate of each user's channel. We formulate 
I each user's QoS requirement as a constraint on the mean square error (MSE) in its received signal, 

■ and we show that these MSE constraints imply constraints on the received signal-to-interference-plus- 

noise-ratio (SINR) of each user Using these MSE constraints, we present a unified design approach 
04 ' for robust linear and non-linear transceivers with QoS requirements. The proposed designs overcome 

T 1 \ 

\ the limitations of existing approaches that only provide conservative solutions and only applicable to 

O ■ 

the case of linear precoding. Furthermore, we provide tractable and computationally-efficient design 

T '— j ' formulations for a quite general model of channel uncertainty that subsumes many uncertainty regions. 

rS ' 

' we also consider the problem of robust counterparts to precoding schemes that maximizes the weakest 
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user's signal subject to a total power constraint on the transmitting antennas. For this problem, we provide 
quasi-convex formulations, for both non-linear and linear transceivers, that can be efficiently solved using 
a one-dimensional bisection search. Our numerical results demonstrate that in the presence of uncertainty 
in the transmitter's knowledge of users' channels, the proposed designs provide guarantees to a larger 
set of QoS requirements than existing approaches, and require less transmission power to satisfy these 
requirements. 
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I. Introduction 

The design of wireless broadcasting schemes that satisfy the quaUty of service (QoS) requirements 
of the intended users (receivers) is of growing interest in interactive communication applications and in 
the downlink of cellular systems with differentiated services. The provision of multiple antennas at the 
transmitter (base station) of the downlink enables the design of schemes that (attempt to) satisfy the users' 
QoS requirements by spatially precoding the users' data in order to mitigate the multiuser interference 
at the (disjoint) receivers. The availability of accurate channel state information (CSI) at the transmitter 
is important in such schemes as it enables the mitigation of the interference experienced by the receivers 
as a result of channel propagation. For scenarios in which one can assume perfect CSI is available at 
the transmitter, the problem of designing a precoder that minimizes the transmitted power required to 
satisfy a set of QoS requirements has been considered in [l]-[6] for the case of linear precoding, and in 
[7]-[ll] for the case of non-linear precoding. 

In practical broadcasting systems, the CSI that is available at the transmitter is subject to a variety of 
sources of imperfection, such as estimation errors, channel quantization errors and short channel coherence 
time. For example, in communication scenarios in which the receivers feed back their quantized CSI to the 
transmitter (e.g., [12]-[14]), the uncertainty in the CSI that is available at the transmitter is dominated by 
quantization errors. Downlink precoder design methods that assume perfect CSI are particularly sensitive 
to these uncertainties, which can result in serious degradation of the quality of the received signals [12], 
[13]. This suggests that the design of downlink precoding schemes should incorporate robustness to 
channel uncertainty. One approach to incorporating robustness is to consider a bounded model for the 
error in the transmitter's estimate of the channels and to constrain the design the precoder so that the 
users' QoS requirements are satisfied for all channels admitted by this model. This bounded uncertainty 
model is useful for systems in which it is difficult to provide the transmitter with an accurate statistical 
model for the channel uncertainty. In particular, it is useful for systems in which users feed back quantized 
channel measurements to the transmitter, as knowledge of the quantization codebooks can be used to 
bound the quantization error. 

For the downlink of cellular systems in which each receiver has a single antenna, the design of a linear 
precoder that minimizes the transmitted power required to guarantee that each user's QoS requirement 
is satisfied for all admissible channels was considered in [15], [16]; see also [3], [17] for designs based 
on a bounded model for the errors in the transmitter's estimate of the (deterministic) autocorrelation 
matrices of the channel. These different approaches approaches formulated the QoS requirements as 
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constraints on the signal-to-interference-plus-noise (SINR) of each user. While the methods proposed 
in [15], [16] provide tractable design formulations and significant improvements in performance over 
previous existing designs, those approaches have two limitations. First, they are not directly applicable 
to non-linear precoding schemes such as Tomlinson-Harashima precoding (THP). Second, when QoS is 
quantified in an SINR sense, the robust linear QoS problem resulted in designs whose tractability is an 
open problem; see also [18]. In order to obtain tractable designs, a conservative design approach was 
taken in [15], [16], and that approach requires the SINR constraints to be satisfied for a superset of the 
original bounded set of admissible channels. In this paper, we provide remedies to both these limitations 
by providing tractable formulations (in the form of semidefinite programs) of both linear and non-linear 
downlink precoding schemes that minimize the transmitted power required to ensure that each user's QoS 
requirement is satisfied for all admissible channels, without expanding the admissible set. We consider 
each user's QoS requirement as a constraint on the mean square error (MSB) in each user's received 
signal, and we show these MSB constraints imply constraints on the received SINR of each user. Since 
the QoS is measured in terms of MSB, our approach is applicable to non-linear Tomlinson-Harashima 
precoding and to linear precoding as a special case. In addition, the proposed designs (for the linear case) 
are obtained with lower computational cost cost than those based on SINR formulations of the QoS in 
[15], [16]. Furthermore, we present a unified treatment of a quite general bounded uncertainty model 
that can represent uncertainty regions resulting from many quantization schemes. The model naturally 
includes channel uncertainty regions that are described using intersection of multiple uncertainty sets, 
e.g., interval constraints on the entries of each user's channel. While we provide exact robust design 
formulations for these types of uncertainties, we also provide conservative formulations that reduce the 
computational complexity of the design for these cases. 

The proposed design approaches can be extended to obtain efficiently-solvable quasi-convex for- 
mulations of some related design problems. In particular, we consider the robust counterpart of the 
problem of maximizing the weakest user's signal (minimizing the largest MSB among the users). For 
precoding schemes that assume perfect CSI at the transmitter, this problem was studied for the case of 
linear precoding schemes in [5], [6]. For the bounded channel uncertainty model, tractable conservative 
approaches to the robust counterpart of the linear minimax precoder design problem were provided in 
[16], but the problem has remained open for the case of non-linear precoding. We provide quasi-convex 
formulations of this robust minimax problem, for both non-linear and linear precoding schemes. These 
formulations can be efficiently solved using a one-dimensional bisection search. We also show that this 
problem can be formulated as generalized eigenvalue problem; e.g., [19]. 
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We also consider the problem of determining the largest uncertainty region for which the QoS 
requirements can be satisfied for all admissible channels using finite transmission power. This problem is 
of considerable interest in the design of quantization codebooks for quantized channel feedback schemes. 
In that case, one might wish to choose the rate of the channel quantization scheme to be large enough (and 
the quantization cells small enough) for it to be possible to design a robust precoder with finite power. We 
provide quasi-convex formulations of this problem, too. Our numerical results demonstrate the efficiency 
of the proposed approaches. In particular, they provide guarantees to a larger set of QoS requirements 
than existing approaches, and require less transmission power in order to satisfy these requirements. 

Our notation is as follows: We will use boldface capital letters to denote matrices, boldface lower case 
letters to denote vectors and medium weight lower case letters to denote individual elements; and 
denote the transpose and the conjugate transpose of the matrix A, respectively. The notation ||x|| 
denotes the Euclidean norm of vector x, while || A|| denotes the spectral norm (maximum singular value) 
of the matrix A, and E{-} denotes the expectation operator. The term tr(A) denotes the trace of matrix 
A, A B denotes the Kronecker product of A and B, and for symmetric matrices A and B, A > B 
denotes the fact that A — B is positive semidefinite. 

II. System Model 

We consider the downlink of a multiuser cellular communication system with Nt antennas at the 
transmitter and K users, each with one receive antenna. We consider downlink systems in which 
Tomlinson-Harashima precoding (THP) is used at the transmitter for multi-user interference pre- 
subtraction. As shown in Fig. 1, the interference pre-subtraction and channel spatial equalization at 
the transmitter can be modelled using a feedback matrix B G ([^KxK ^ feedforward precoding matrix 
P G C^'^^. Since linear precoding is the special case of the THP model in which B = 0, we will 
focus our development on the THP case and will extract the special case results for linear precoding as 
they are needed. The vector s G in Fig. 1 contains the data symbol destined for each user, and we 
assume that is chosen from a square QAM constellation S with cardinality M and that E{ss^} = I. 
The Voronoi region of the constellation V is a square whose side length is D. 

In absence of the modulo operation, the output symbols of the feedback loop in Fig. I, Vk, would be 
generated successively according to = Sk — Z]^=i ^k,jVj, where at the A;* step, only the previously 
precoded symbols vi, ..,Vk-i are subtracted. Hence, B is a strictly lower triangular matrix. The role of 
the transmitter's modulo operation is to ensure that Vk remains within the boundaries of V, and its effect 
is equivalent to the addition of the complex quantity = i^^f D + j D to v^, where i^^ , j^^'^d g 
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Fig. 1. Multiple-input single-output downlink system with Tomlinson-Harashima precoding at the transmitter. 
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Fig. 2. Equivalent linear model for the transmitter. 



and j = V— 1 ■ Using this observation, we obtain the standard linearized model of the transmitter that 
does not involve a modulo operation, as shown in Fig. 2; e.g., [20]. For that model, 

v = (I + B)-iu, (1) 

where u = i + s is the modified data symbol. As a result of the modulo operation, the elements of v 
are almost uncorrelated and uniformly distributed over the Voronoi region V [20, Th. 3.1]. Therefore, 
the symbols of v will have slightly higher average energy than the input symbols s. (This slight increase 
in the average energy is termed precoding loss [20].) For example, for square M-ary QAM we have 
a^, = E{|vfcp} = j^j^E{|sfcp} for all k except the first one [20]. For moderate to large values of M 
this power increase can be neglected and E{vv^} = I is often used; e.g., [20]-[22]. Hence, the average 
transmitted power constraint can be written as Ev{x^x} = tr(P^P). 
The signals received at each user, y^, can be written as 

yk = hfcx + nfc = hfcP(I + B)"^u + nfc, (2) 
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where h/t G C^^^* is a row vector representing the channel gains from the transmitting antennas to the 
A;* receiver, and Uk represents the zero-mean additive white noise at the A;* receiver, whose variance is 
0"^^. At each receiver, the equalizing gain is used to obtain an estimate = (7fchfcP(I + B)^^u + (7fcnfc 
of the modified data symbol u^. Following this linear receive processing step, the modulo operation is 
used to obtain s by eliminating the effect of the periodic extension of the constellation caused by the 
integer ik- In terms of the modified data symbols, we can define the error signal 

Uk-Uk = (^fchfcP - rrifc - bfc)v + g^Uk, (3) 

where rrifc and bfc are the A;* rows of the matrices I and B, respectively. The error signal in (3) is 
equivalent to — when the integer is perfectly eliminated. This occurs with high probability even 
at reasonably low SINRs. Using this error signal, the Mean Square Error (MSE) of the /c* user is given 
by 

MSEfc = E{|?ifc-Ufcp}= llgfchfcP - mfc - bfcf + Ifffcpcr^ 



2 

5fehfcP - rrifc - bfc g^ak 



III. Transceiver Design with MSE Constraints: Perfect CSI case 



(4) 



We consider downlink scenarios in which each user has a quality of service constraint on the form of 
an upper bound on its mean square error, MSE^. The formulation of QoS design problem in terms of 
the MSEs is motivated by the following result. 

Lemma 1: For any given set of uses' channels h^, if there exists a transceiver design P,'B, gk that 
guarantees that MSE^ < Cfc> then that design also guarantees that SINR^ > (1/Cfc) — 1- 

Proof: See Appendix A. ■ 
It worth observing that Lemma 1 provides a relation between the guaranteed MSE of each user and 
its achievable SINR that extends to scenarios in which users' channels are not accurately available at 
the transmitter. This is valid as long as long as the MSE is guaranteed for all channels within a given 
uncertainty set around the transmitter's estimate of the users' channels. Analogous relation existed for 
when accurate CSI is assumed to be available at the transmitter, e.g., [9]. In order to facilitate our 
development of robust precoding schemes with QoS constraints, we will briefly consider the design 
problem for the case of accurate transmitter's knowledge of the users' channels. In that case, the design 
of the downlink transceiver components P, B and gk so as to minimize the total transmitted power 
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subject to satisfying the users' MSE requirements can be formulated as 



mm 
subject to 



|vec(P) 



(5a) 
(5b) 



Since the norm in (5b) is unitarily invariant, problem is independent on the phase of each g^. Indeed, 
if {\gk\ e^^''} and P are the optimal equalization gains and precoding matrix, respectively, then 
and P Diag(e-'^S . . . , e-'^^ ) are also optimal. Hence, there is no loss of generality in choosing the all the 
equalization gains to be real. Using this observation and the definitions 



hk 



mk 

fk 



Re{hk} Im{hfc} 

Re{P} Im{P} 
- Im{P} Re{P} 

Re{hk}/gk lm{hk}/gk 
Re{mfc} Im{mfc} 



(6) 

(7) 

(8) 
(9) 
(10) 



the design problem in (5) can be formulated as a convex Second Order Cone Program (SOCP) 



min t 

P, B, fk,t 



subject to ||vec(P)|| < t, 

II [hfcP - fkmk - bfe, o-„J II < \/Cfc/fc l<k<K, 



(11a) 

(lib) 
(11c) 



When the channel is accurately known at the transmitter, we can obtain stronger conclusions than those 
provided by Lemma 1. In this case, it can show be shown using a contradiction argument that solution 
to (11) results in MSEfc = Cfe> the constraints in (11c) are achieved with equality at optimality. Since 
at optimahty, each MSE^ = 1/Cfc, we have that for the perfect CSI case, SINR^ = l/MSE^ - 1, e.g., 
[9]. 

Another advantage of the convex conic formulation in (11) is the possibility to include shaping 
constraints (e.g., [23]) on the transmitting antennas. These constraints are expressed as either second 
order cone or positive semidefiniteness constraints on the precoding matrix P. More importantly, the 
convex formulation in (11) enables us to derive robust counterparts of the original design problem in (5) 
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for uncertainty models presented in the following section. 

IV. Channel Uncertainty Model 

We will consider an additive uncertainty model on the form: 

J 

Uk{h) = {hfc I hfc = hfc + efc = hfc + Wj^j, w^Qw < ,5^}, (12) 

where is the transmitter's estimate of the A;* user's channel, and is the corresponding estimation 
error. The above model enables us to treat several different uncertainty regions in a unified way. For 
example, it can model the following uncertainty sets: 

• Ellipsoidal and Spherical Uncertainty Sets: By choosing Q = I, the uncertainty set in (12) describes 
an ellipsoidal uncertainty region around the channel estimate h^. The spherical uncertainty set with 
center h;. and radius 5k is the special case that arises when cf)j are selected to be the rows of l2Nt 

• Interval Uncertainty Sets: Interval constraints on each element of can also be modeled as 
uncertainty sets of the form in (12). By taking cf)j to be the rows of l2Nt ^^d Q to be the matrix 
whose only non-zero element is Qa = 1, then the uncertainty set in (12) models an interval constraint 
on the entry of the error h;^.. Interval constraints on multiple entries of h;;. can be represented as 
the intersection of uncertainty sets on the form (12); see Section V-A. 

This additive uncertainty model is useful for systems in which the channel state information is quantized 
at the receivers and fed back to the transmitter; e.g., [12]. If a vector quantizer is employed at the 
receivers, then the quantization cells in the interior of the quantization region can be often approximated 
by ellipsoids. On the other hand, if a simple scalar quantizer is employed, the quantization regions can 
be modeled using a set of interval constraints. 

V. Transceiver Design with MSE Constraints: Uncertain CSI Case 

In this section, we will design a robust transceiver that minimizes the total transmitted power necessary 
to guarantee that the users' MSE requirements are satisfied for all admissible channels h^. in the 
uncertainty region Uk{6k) in (12). Using the formulation in (11), this robust problem can be stated 
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as 

min t (13a) 

P,B, U, t 

s. t. ||vec(P)|| < t, (13b) 
II [hfcP - Amfc - bfc, a„j|| < Ta/fc V hfc G 1 < A: < JC. (13c) 

This is a semi-infinite conic programming problem. In particular, the constraint (13c) represents K infinite 
sets of second order cone (SOC) constraints, one for each hj^. G Uk{5k). However, we can precisely 
characterize each of these infinite sets of SOC constraints using a single Linear Matrix Inequality (LMI), 
as stated by the following theorem. 

Design Formulation 1: The robust transceiver design problem in (13) is equivalent to the following 
semidefinite program (SDP) 



min t 
fj,. t 

s. t. ||vec(P)|| < t, 



Ckfk - fJ-k 











[hfcP - mfc/fc -bfc, an. 



AifcQfc 4[P, 0] 



l<k<K, 



(14a) 
(14b) 

>0, 
(14c) 



□ 

Proof: See Appendix B. ■ 
This result shows that the original design problem in (13) with an infinite set of constraints is equivalent 
to the convex SDP in (14) that can be efficiently solved using interior point methods, e.g., [24]. Such 
equivalence is an advantage of the structure of the uncertain parameter of the SOC representation, in 
(13c). In these SOC constraints, the channels h^., and consequently the uncertain parameters only exist 
on one side of the SOC. Hence, exact characterization of these SOC with uncertain parameters can be 
obtained. In contrast, when the QoS requirements are of the form of bounds on the SINR, then even in the 
case of linear precoding, both sides of the SOC constraints that enforce the QoS requirement depend on 
h^, and the resulting design problem is not known to be tractable [18, pp. 7]. In [15], [16] this unknown 
tractability was addressed by taking a conservative approach to the robust design problem. By adopting 
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MSE constraints, the QoS robust design problem can be efficiently solved. 

A. Case of Intersecting Uncertainty Sets for each 

The design problem in (13) and its efficiently-solvable equivalent in (14) were obtained for an 
uncertainty set of the form in (12). However, that design extends naturally to the case in which the 
uncertainty region for each is described as the intersection of more than one uncertainty set of 
the form (12); that is, uncertainty set is of the form 

L 



Uk = f] Ui{5k) (15) 



Note that there is no restriction in assuming that each has the same uncertainty size parameter 6k, 
since Q^, in (12) can be chosen to accommodate different sizes and geometrical regions. Examples of 
constraint sets of the form in (15) include interval constraints on each of the entries of discussed in 
Section IV. In that case, there is a constraint on the form of (12) for each element of h^,, and hence 
L = 2Nt. To formulate the modified version of (13) with replaced by as an SDP, one simply 
observes that each Ul{5k) in (15) constitutes an LMI of the form in (14c) and that all LMIs must be 
satisfied. Hence, the resulting SDP has LK LMIs. While that formulation is precise, when L is large it 
might be prudent to adopt a conservative formulation with fewer constraints. This conservative approach 
involves enveloping (15) in a superset that can be described more efficiently, and then requiring the MSE 
constraints to be satisfied for all channels in this superset. The following theorem describes a conservative 
approach that results in a formulation with the same number of LMIs as that in (14). This represents a 
reduction in the number of LMIs by factor of L. 

Design Formulation 2: The solution of robust transceiver design problem in (13) for the intersection 
of uncertainty sets in (15) is upper-bounded by the solution of the following SDP 

min t (16a) 
s. t. ||vec(P)|| < t, (16b) 



[hfcP - fkmk - hk, <^nf of Vafki 



> 0, 



l<k<K. (16c) 
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□ 

Proof: See Appendix B. ■ 

B. Largest Feasible Uncertainty Size 

In this section we consider the related design problem of finding the largest value of the uncertainty 
size 5, namely 5max for which there exists a robust transceiver of finite power that satisfies the MSE 
constraints for all admissible channels in the uncertainty region of size 5max- This problem is connected 
to the problem of designing codebooks for the quantization of the users' channels. The codebook design 
needs to yield quantization regions that can be "covered" by uncertainty sets of size 5max in order for 
the robust transceiver design problem to be feasible. Using the problem formulation in (14), finding the 
value of 5max is equivalent to solving 

max p (17a) 

P, B, A-,M,P 

s. t. Afc(a,/9)>0, l<k<K, (17b) 

where Afc(Cfc, p) is defined in (14c). Since p is an optimization variable rather than a design parameter, 
the bilinear terms in Ak{(k, p) mean that design problem is not jointly convex in the design variables p 
and P. However, the problem is quasi-convex (c.f. [25]), and an optimal solution can be efficiently found 
using a one-dimensional bisection search on p in which the problem solved at each step is the convex 
feasibility problem corresponding to (17) with a fixed value for p. For the case of the intersection of 
uncertainty regions in (15) the constraint ^kiCk^p) in (16c) can be used in the place of (17b). In that 
case, the optimal value of the design problem becomes a lower bound on 6max- It is worth observing that 
largest uncertainty size for the special case of linear precoding is less than of its THP counterpart. This 
follows by observing that finding 5max in the linear precoding case solves a restriction of the problem 
(17) in which each is set to 0. This result does not use the assumption of negligible precoding loss, 
since 6m.ax is the maximum uncertainty supported by any transceiver design with finite power. 

VI. Robust Counterpart of Fair Minimax Transceiver Design 

In the previous section, the focus was on the robust counterpart of the transceiver design problem that 
minimize the total transmitted power subject to satisfying the users' MSE constraints. In this section, we 
consider the related problem of minimizing the maximum MSE among aU users subject to a transmitted 
power constraint, in scenarios with uncertain CSI. This design problem provides a notion of fairness 
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among users based on the value of their MSEs. While this problem has been considered in scenarios that 
assume perfect CSI in [5], [6], we can formulate the robust counterpart of this design problem under the 
channel uncertainty model in (12) as the following semi-infinite quasi-convex optimization problem 



mm 



(18a) 



s-t. ||[hfcP-/fcmfc-bfc, a„j||<7Co/fc, Vh;,EWfc(4), I < k < K, (18b) 



-tr(PP^) < Ptotal. 



(18c) 



Using the characterization in (14c) of the infinite set of SOC constraints in (18b), this design problem 
can be formulated as the following quasi-convex optimization problem 



mm 



s. t. Afc(Co,<^fc) > 0, l<k<K, 
||vec(P)|| < V2 Ptotal. 



(19a) 

(19b) 
(19c) 



This problem can be efficiently solved by using a bisection search on in which problem solved at 
each step is the convex feasibility problem generated by (19) with a fixed value of Furthermore, 
we can observe that each constraint in (19b) can be written as 



i/gk 

- iii,Ou SulP. 01 > 0. 

1/gkI 



VCkfk - f^k [hfcP - fkmk - hk, (^n,. 

/XfcQfc 6k[P, 0] 

[hfcP - /fcmfc - bfc, cr„J (^fe[P, 0] VCfc/fcl 

Hence, (18) is equivalent to minimizing the largest generalized eigenvalue of a pair of (block diagonal) 
symmetric matrices that depend affinely on the decision variables [19], [26] — a problem that takes the 
form 



mm a 

X,CJ 



s. t. qA^(x) - A2(x) > 0, 
Ai(x) > 0. 
B(x) > 0. 



(20a) 

(20b) 
(20c) 
(20d) 
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this observation allows us to employ efficient algorithms that exploit the structure of the constituent 
matrices in (??) are available for such problems; c.f. [19], [27]. 

VII. Simulation Studies 

In this section, we demonstrate the performance of the proposed robust QoS designs for non-linear 
Tomlinson-Harashima precoding (RTHP-order 1, 2) and linear precoding (RLin) that were presented in 
Section V. For Tomlinson-Harashima precoding, ordering of the users' channels is necessary prior to 
precoding. Finding the optimal ordering requires exhaustive search over all possible permutations of the 
transmitter's estimate of users' channels h;., and instead of that we have implemented two suboptimal 
ordering methods. The first method applies the BLAST ordering in [28] to the transmitter's estimate of 
users' channels. The second method is a generalization of the ordering method in [29] that selects a 
channel ordering that minimizes the reciprocals of the received SINRs when the precoder matrix P is 
an identity matrix. In our generalization, the ordering selection criteria is minimizing the sum of each 
user's SINR requirements divided by its received SINR (when P = I), a quantity that is proportional to 
the power necessary for each user to achieve its SINR requirement. 

In our numerical studies, we consider a spherical uncertainty region of radius 6k for each user. This 
model will facilitate the comparisons with other existing approaches for the linear precoding model, 
namely the robust autocorrelation matrix approach in [3], [4] (Robust Correl. Appr.), the robust power 
loading approach (RLin-PLl) using SINR constraints in [17], and the robust power loading approach 
(RLin-PL2) using MSE constraints in [30]. We will also compare with the conservative approach to robust 
linear precoding with SINR constraints in [15], [16]. The work in [15], [16] presented three conservative 
approaches and we are comparing with the best conservative approach, namely the "Structured SDP" 
approach in Section VI of [15]; see also Section IV [16]. As we make the comparisons, we would like 
to point out that all previous approaches do not extend to non-linear Tomlinson-Harashima precoding, 
but the approaches proposed herein are inherently applicable to both linear and Tomlinson-Harashima 
approaches. In order to totally specify the schemes used in our comparisons, we point out that the 
approaches in [17] and [] requires the beamforming vectors (normalized columns of P) to be specified. 
We will use the zero-forcing beamforming vectors (the columns of the pseudo-inverse of H). In addition, 
the approaches in [3], [4] and [17] are based on uncertainty models that are different from the one in 
(12), and from each other. The approach in [3], [4] considers a model in which the spectral norm of 
the error in the (deterministic) autocorrelation matrix = h^h^ is bounded, and in the approach 
in [17] the Frobenius norm of the error in is bounded. However, by bounding these norms of 
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in terms of the norm of e^, a comparable uncertainty set can be generated. We will compare these 
schemes in an environment with Nt = 3 transmit antennas and K = 3 users. In our experiments, we will 
evaluate performance statistics for the standard case of independent Rayleigh fading channels in which the 
coefficients of the fading channels are modeled as being independent circular complex Gaussian random 
variables with zero-mean and unit variance, and the receivers' noise sources are modeled as zero-mean, 
additive, white, and circular Gaussians with unit variance. 

A. Performance Comparisons against SINK Requirements 

In this comparison, we randomly generated 2000 realizations of the set of channel estimates {h^}^]^ 
and examined the performance of each method in the presence of uncertainties of equal sizes, 6k = 6 = 
0.05, VA;. The SINR requirements of the three users are also equal. For each set of channel estimates and 
for each value of the required SINR we determined whether each design method is able to generate a 
precoder (of finite power) that guarantees the required SINRs. In Fig. 3 we plot the fraction of the 2000 
channel realizations for which each method generated a precoder with finite power against the equal 
SINR requirement of each user. 

From this figure, it clear that the proposed robust designs for linear (RLin) and non-linear (RTHP- 
order 1, 2) precoding satisfy the SINR requirements for larger percentages of channels. The robust 
conservative approach for linear precoding (RLin-Conservative) [15], [16] and the power loading method 
in [] achieve the QoS requirements for a percentage of channels that is quite close to that of the proposed 
linear approach (RLin). For the robust linear power loading approach (RLin-PL2) in [30], the QoS design 
problem in terms of MSB constraints was justified as a heuristic measure for the SINR requirements. 
However, using Lemma 1 we showed that the MSB constraint of each user implies a minimum achieved 
SINR. 

For the comparison in Fig. 4, we selected all the sets of channel estimates from the 2000 sets used 
in the previous experiment for which all design methods were able to provide robust QoS guarantees 
for all SINRs less than or equal to 6dB, and we calculated the average, over the 274 such channel 
environments, of the transmitted power required to achieve these robust QoS guarantees. We have plotted 
the equal SINR requirement of each user versus the average transmitted power in Fig. 4. The average 
transmitted power approaches infinity for a certain value of SINR when for one (or more) of the channel 

'a bound on the spectral norm of the error in the matrix can be obtained as follows: ||(hfc + efc)^(hfc +efc) — h^hfc|| — 
Wh^Bk + hfe + e^efcll < |lhf efe|l + ||ef hfe|| + ||efefc|| = 2||hfc |1 ||efc |[ + ||efef . The same bound also holds for the 
Frobenius norm, since the matrices on the immediate right hand side of the inequahty are all rank one. Furthermore, the 
uncertainty — 5fchfc/||hfc|| achieves this upper bound with equality for both norms. (See also [31].) 
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Fig. 3. Percentage of the 2000 channel realizations for which the robust QoS guarantee can be made against the required 
SINRs, for a system with Nt = 3 and K = 3. 



estimates the method under consideration cannot provide the robust QoS guarantee with finite power. 
This figure demonstrates the saturation effect that channel uncertainty imposes on the growth of the 
SINR of each user with the transmitted power for both of Unear and non-Unear precoding. This effect 
was observed in [12] for non-robust linear precoding on the MISO downlink with quantized CSI. Fig. 4 
also illustrates the role that robust precoding can play in extending the SINR interval over which linear 
growth with the transmitted power can be achieved. This is particularly evident for the robust non-linear 
approaches (RTHP-order 1, 2) and the robust linear approach (RLin). We also observe that the second 
ordering method for Tomlinson-Harashima precoding provides better performance than the first one, since 
it selects the channel ordering in a way that attempts to minimize the sum of powers necessary to achieve 
each SINR requirement. 

B. Performance Comparisons against Uncertainty Size 

In this comparison, we used the 2000 randomly generated realizations of the set of channel estimates 
{^k]k=i examine the performance of each method in the presence of equal uncertainty, 5^ = 5, 
VA;. The QoS requirement of each user is such that the SINR is at least 10 dB. In Fig. 5 we provide 
the percentage of the 2000 channel realizations for which each method generated a precoder with finite 
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power as a function of the size of the uncertainty. From this figure, it is clear that for a large range 
of uncertainty sizes, the proposed non-linear approaches (RTHP-Order 1, 2) provide SINR requirements 
for many more channel realizations than other approaches. This is due to the fact that the proposed 
linear approach is a special case of the proposed THP design, and the other existing linear approaches 
are either conservative or restricted to the optimization of powers for a given transmit directions. While 
the conservative linear precoding approach (RLin-conservative) in [15], [16] and the robust linear power 
loading approaches (RLin-PL2) in [30] have close performance in terms of number of channel realizations 
for which the methods achieves the QoS requirements, they do use more power in order to achieve the 
QoS requirements as shown in Fig 6. 

In Fig 6, we selected those sets of channel estimates from the 2000 sets used in the previous experiment 
for which all design methods were able to provide robust QoS guarantees for all uncertainties with 
6 < 0.015. We calculated the average, over the 614 such channel environments, of the transmitted power 
required to achieve these robust QoS guarantees and we have plotted the results for different values of 
6 in Fig. 6. The average transmitted power approaches infinity for a certain value of 6 when for one (or 
more) of the channel estimates the method under consideration cannot provide the robust QoS guarantee 
with finite power. It is clear from Fig. 6 that the proposed robust Tomlinson-Harashima designs are 
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Fig. 5. Percentage of channel realizations for which the robust QoS guarantee can be made against the uncertainty size S, for 
a system with Nt — 3 and K = 3. 



capable of achieving SINR requirements for larger values of uncertainty sizes than other approaches. 



We presented a unified design approach for robust non-linear and linear transceivers with users' QoS 
requirements subject to deterministically-bounded channel uncertainty model. The proposed approach 
formulated the QoS requirements in terms of MSE constraints and showed that these constraints imply 
corresponding constraints on the achieved SINR of each user. It provided (convex) semidefinite program 
formulations of the design problem that are efficiently-solvable. Furthermore, these design formulations 
were obtained for a quite general model of channel uncertainty that include many uncertainty regions. We 
also showed how these designs can be used to provide quasi-convex designs for the robust counterpart 
of the problem of fair transceivers that maximizes the signal quality of the user with the weakest signal. 
Numerical results demonstrated that under uncertain CSI conditions, the proposed designs provided 
guarantees to a larger set of QoS requirements than existing approaches, and require less transmission 
power to satisfy these requirements. 



VIII. Conclusion 
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Appendix A 
Proof of Lemma 1 

Consider the quantity (uk — Uk) in equation (3). Assuming correct removal of i^, we have 

Sk- Sk = (s'fchfcP - Hifc - bfc)v + gkUk, 

or equivalently, 

Sk = akSk + ^ a-iSi + aoUk, 
where 1^ is the set of interfering symbols with s^. Using (22), we can write 

ieik 

E{|sfc-SfcP} = kfc - + ^ lojp + laolcTfc, 

= E{|4|'} + l-2Re{afc}, 
l + l/SINRfc = E{|sfc|2}/|afc|2. 
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Consider the MSE constraint E{\sk — SfcP = E{|sfcp} + 1 — 2Re{afc}} < Ck < ^- This can can be 
written as 

E{|sfcp}(l-Cfc) < 2Re{afc}(l-a)-(l-a)' (27) 

= Re2{afc}- (Re{afc}}-(1-Cfc))' < lafcp. (28) 

The latter inequality is equivalent to 1 + 1/SINRfc < 1/(1 - Cfc)> or equivalently SINRfc > (1/Cfc) - 1- 

Appendix B 
Proofs of Design Formulations 1 and 2 

The proofs are based the following lemma which is a concatenation of two results in [32]: 
Lemma 2: Consider the SOC constraint ||Ax + b|| < y for every [A, b] in the uncertainty region 
given by 

J 



U 



[A, b] I [A, b] = [A", bO] + % [A^ b-'l, eev] 
I 0^ 0, / = 1,... (29) 



Then the set Si of pairs (x, y) satisfying || Ax + b|| < y for every [A, b] e Z// is subset of the set ^2 
of pairs (x, y) such that there exist non-negative scalars fi^, . . . ,fi^ satisfying 



y-Ei/ (AOx + bO)^ 

Ei/Q^ [A^x + bi ... A^x + b^]^ 

A^x + b [A^x + b^ ... A'^x + b'^] yl 



> 0. (30) 



When L = 1, 5i = 52- □ 
To prove Design Formulation 1, we involve Lemma 2 with L = 1 to show the equivalence between 
the SOC constraints in (13c) and the corresponding LMIs in (14c). The nonnegativity constraints on each 
fik is implied by positive semidefiniteness of the diagonal blocks of the matrices in (13c) The proof of 
the Design Formulation 2 is similar, but when L > 2 the application of Lemma 2 results in a conservative 
design formulation, and hence an upper bound on the required transmission power. 
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